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PACS 82. 70. -y - Disperse systems; complex fluids. 

PACS 61.20.Qg - Structure of associated liquids: electrolytes, molten salts, etc. 

Abstract. - We predict the nature (attractive or repulsive) and range (exponentially screened or 
long-range power law) of the electrostatic interactions of oppositely charged and planar plates as 
a function of the salt concentration and surface charge densities (whose absolute magnitudes are 
not necessarily equal) . An analytical expression for the crossover between attractive and repulsive 
pressure is obtained as a function of the salt concentration. This condition reduces to the high-salt 
limit of Parsegian and Gingell where the interaction is exponentially screened and to the zero salt 
limit of Lau and Pincus in which the important length scales are the inter-plate separation and 
the Gouy-Chapman length. In the regime of low salt and high surface charges we predict - for 
any ratio of the charges on the surfaces — that the attractive pressure is long- ranged as a function 
of the spacing. The attractive pressure is related to the decrease in counter-ion concentration as 
the inter-plate distance is decreased. Our theory predicts several scaling regimes with different 
scaling expressions for the pressure as function of salinity and surface charge densities. The 
pressure predictions can be related to surface force experiments of oppositely charged surfaces 
that are prepared by coating one of the mica surfaces with an oppositely charged polyelectrolyte. 
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Introduction. — The interactions between oppositely 
charged surfaces are important in both biological and ma- 
l^erials science contexts. Proteins contain both cationic 
and anionic regions and in some cases, interactions be- 
liween proteins are due to unlike charge attraction, medi- 
ated by the intervening counterions and salt. The de- 
livery of cationic vesicles to cells - of interest in gene 
therapy applications - involves counterion and salt me- 
diated electrostatic interactions of two oppositely charged 
membranes [1~3]. Similar considerations may also be im- 
portant in understanding adhesion and fusion in systems 
of oppositely charged bilayers [4,5]. Recent experiments 
on hydrophobically prepared mica surfaces have indicated 
that such surfaces have domains with different charges 
and the observed long-range attractions (in the nanometer 
regime) may again be related to unlike charge attractions 
mediated by counterions and salt [6,7]. 

In this paper we predict the interactions between two 
homogeneously charged surfaces with opposite charge. 
The surfaces are in aqueous solution that contains the 
counterions and added salt. The interactions can be at- 
tractive or repulsive, short-ranged (exponentially decay- 
ing) or long-ranged (power-law) depending on the ratio 



of the distance between the surfaces to the important 
length scales of the problem: (i) The Gouy-Chapman 
length, Agc = l/(27rZB(T) that is inversely proportional 
to the surface charge density cr; (ii) the Bjerrum length 
Ib = e'^/eftBT equal to about 7 A in water (e ~ 80) 
at room temperature; and, (iii) the Debye-Hiickel (DH) 
length. Ad — l/V^TrlsCb where ct is the bulk 1:1 salt con- 
centration. Since there are several length scales there are 
several regimes that characterize the interactions. 

In order to model interactions between charged surfaces, 
membranes and particles, one typically considers two pla- 
nar surfaces separated by a distance, d. Previous studies 
considered the symmetric case where the two surfaces have 
fixed and equal surface charge densities. Within Poisson- 
Boltzmann (PB) theory, it can be shown that the interac- 
tion between two surfaces with the same charge is always 
repulsive due to the counter-ion entropy. Other refine- 
ments include correction to the PB theory especially in the 
limit of strong surface charges and multi-valent counteri- 
ons. The theory of the interactions between two surfaces 
with opposite charges has received less attention. The pi- 
oneering study of Parsegian and Gingell [14] considered 
high salt concentrations and used the linear DH theory, to 
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predict regions of repulsive and attractive inter-plate pres- 
sure as a function of the distance and salt concentration. 
Lau and Pincus [13] considered the case of two oppositely 
charged surfaces with no added salt and found a simple 
analytical criterion for the crossover between attraction 
and repulsion. 

The interesting, intermediate region of high surface 
charge and relatively low salt concentrations requires the 
use of the non-linear PB theory. There can be a large 
regime in which the distance between the plates is smaller 
than the DH screening length (which can be very large 
at low salt concentrations) but larger than the Gouy- 
Chapman length (which can be very small for high sur- 
face charge) in which the non-linear treatment must be 
applied. The case of two plates with equal and oppo- 
site surface charges was considered recently by several au- 
thors [15,16] . In this case, where the two plates are electri- 
cally neutral and the ions in solution have no net charge, 
the force is always attractive. At some characteristic inter- 
plate distance, d*, which scales as Ago ln(AD/AGc), the 
counterions arc released into the; bulk reservoir provided 
that Ago d* Ad- In addition, for inter-plate sep- 
aration in the range of intermediate values (between Ad 
and Ago) I the pressure was shown to be long-ranged, and 
scales with the inverse of (P [15]. In the present work, 
we generalize these results and consider two oppositely 
charged surfaces whose charges are not necessarily equal 
in magnitude. This general case is important in order to 
analyze experiments in which surfaces are not completely 
antisymmetric. We use our analytical theory to predict a 
crossover between attraction and repulsion and the coun- 
terion release concept is extended to the asymmetric case. 
We also are able to consider several limiting regimes for 
the asymmetric case that complement the numerical solu- 
tions of the problem. 

Poisson Boltzmann model. — The Poisson- 

Boltzmann (PB) theory is a mean-field theory that relates 
the electric potential, V'(0> ^'^d the Boltzmann distribu- 
tion for the ion number density, c(r), at thermodynamic 
equilibrium. For two surfaces immersed in a 1:1 mono- 
valent ionic solution, and for a dimensionless potential 
(f) = eil)/k-B.T, the PB equation reads: 

= Kq sinh0 = Aq^ sinh(/) , (1) 

where Ad = ^d^ is the Debye-Hiickel (DH) screening 
length defined above. 

We consider two charged surfaces that are infinite in 
extent in the [x, y) plane and are separated in the z di- 
rection by a distance d. In this case, the PB equation 
reduces to a one-dimensional equation in the coordinate z 
perpendicular to the planes. The focus of this paper is on 
two oppositely charged plates that we call the asymmetric 
two plate problem. A positively charged plate with charge 
density (T+ > is located at z = d/2, and a negatively 
charged one with (t_ < is located at = —d/2. While 
in the well-studied symmetric case, where C7+ = C7_, [8] 



the potential (and ion densities) are symmetric about the 
midplane ^; = 0, in the asymmetric case, the midplane is 
no longer a plane of symmetry. Instead, a separate bound- 
ary condition at each plate must be explicitly considered. 
These conditions relate the electric field at each plate to 
the surface charge density: 

^ =47r/B|a±|^^ >0 (2) 

z=±d/2 -^i 

where A± are the Gouy-Chapman lengths for the corre- 
sponding surfaces. 

The spatial dependence of the potential and ion densi- 
ties is obtained by solving the PB equation, eq. (1) subject 
to the boundary conditions, cq. (2). The profiles predict 
the local concentration of the mobile ions and their associ- 
ated potentials; these can be measured, for example, using 
scattering techniques. However, more often, the forces ex- 
erted on the charged plates are measured [9]. For a given 
separation d, the pressure (or equivalently, the force per 
unit area) must be constant in the entire region between 
the plates if the system is in thermodynamic equilibrium. 
Thus, it is only necessary to calculate the pressure at any 
convenient point, z, within the gap. Because the two-plate 
system is in contact with a reservoir of mobile ions (the 
salt reservoir), the net pressure (in units of U-qT) exerted 
on the plates is given by the difference between the inner 
and outer pressures, H = Pin — fout- This pressure can 
be calculated, for example, by integrating the PB equa- 
tion once and relating the integration constant with the 
pressure H [8]: 

The z-independent osmotic pressure comprises two terms: 
(i) an attractive contribution whose origin is the electro- 
static energy; this has the form of a negative term pro- 
portional to the square of the electric field [10] and (ii) 
a repulsive contribution that arises from the translational 
entropy of the ions and is given by the ideal-gas law. In 
this term we have already subtracted the outer pressure, 
-Pout = 2c(, as is explained above. 

It is easy to show that the pressure for the case of 
symmetrically charged plates is always repulsive (H > 0) 
within the PB approximation. However, in the gen- 
eral, asymmetric case the pressure can be either repulsive 
(H > 0) or attractive (H < 0) as we explain below. 

Our two-plate problem is fully determined by four phys- 
ical parameters: The two surface charge densities a± , the 
ionic strength, c;,, and the separation d. However, by using 
normalized variables, it is easy to show that the problem is 
uniquely defined by three ratios: A±/Ad and d/Xu, where 
A± is related to i7± in eq. (2) and Ad is related to ci,. In our 
PB solution we use an alternative parametrization scheme 
including the dimensionless pressure H = H/cft, the sur- 
face potentials (f)± = (j){±d/2). We can now relate those 
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three parameters to the original ones by three relations. 
Two relations can be derived from eq. (3): 



n = 2(cosh(/.±-l)- ( ^ 



(4) 



Integration between the two boundaries, ±d/2, gives an 
additional, third relation: 



_d_ '■'t>+ 
Ad 



J A. 



2(cosh0'-l)-n 



(5) 



We can now express the PB solution and pressure 11 [via 
eqs. (4) and (5)] as a function of the parameters that char- 
acterize the physical system namely, the surface charge 
densities, A_|-,A_, the DH length, Ad (determined by the 
salt concentration), and the separation between the plates, 
d. For a few simplified cases, an analytical solution exists, 
while in the general asymmetric case wc can write the 
potential, 0, in terms of elliptic integrals [11] whose solu- 
tion can be obtained only numerically. We next present 
a general analytical result that predicts when the inter- 
action crosses over from repulsive to attractive as a func- 
tion of the system parameters. The physical origin of this 
crossover is the competition between the electrostatic and 
entropic interactions described above. 

Attractive to repulsive crossover. — The condi- 
tion n = in eq. (3) for the asymmetric, two plate system 
determines the cross-over from repulsive to attractive in- 
teractions in the system: 



(<^') +2c6(cosh<?i- 1) = 0. 



(6) 



This is a relation between the potential. (p{z), and its 
derivative <j)' at any point z. The condition, 11 = 0, also 
fixes one relation between the three dimensionless ratios: 
A±/Ad and d/Xu- Namely, this confines the system to a 
two dimensional surface in the three-dimensional parame- 
ter space (A+/Ad, A_/Ad, ^/Ad). 

An analytical expression for this crossover is found by 
observing that the case of 11 = can be c;xactly mapped 
onto the equations that describe a single plate in contact 
with the same reservoir [12]. For this purpose, we con- 
sider a system with a single, positively charged plate at 
z — 0. The analytical expression for the potential and its 
derivative (electric field) are well known 
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where 7+ = \J (A+/Ad)^ + 1 - (A+/Ad). The mapping 
between the two problems is simply done by requiring that 
at distance d away from the z = charged plate (the single 



plate case) the electric field is equal to the electric field as 
determined from Gauss's law at the negative plate located 
at ^; = —d/2 in the two-plate problem: (j)'{d) = — 2/A_. 
This results in the relation: 



1+ 



(9) 



where 7_ is similarly defined as 7_ = ij (A_/Ad)^ + 1 — 
(A_/Ad). 

The relation (9) between 7-j- (or X±) is equivalent to 
the n = crossover in the asymmetric two plate system, 
between the repulsive, 11 > 0, and attractive, 11 < 0, 
regimes. When 7_ ^ (the negative plate is neutral, 
c7_ 0), the plates must repel each other. In addition, 
although we have only treated so far the case it+ > \a-\, 
our results are quite general since the two-plate pressure is 
invariant under the exchange a+ )fT_). Therefore, the 
condition for attraction reads: 



7- 



(10) 



This result is plotted on fig. 1 where two lines separate a 
central region of attractive interactions from two wedges 
in the (<j-, (t+) plane, that denote repulsive interactions. 
In our plots, the charge densities are normalized hy ao = 

l/2lTl-Bd. 

This result is exact for arbitrary salt concentration and 
surface charge densities. It has two limits that have been 
previously studied. One limit, that of zero salt, was an- 
alyzed by Lau and Pincus [13]. In this limit, the coun- 
terions in the solution balance the surface excess charge 
Act = (T+ — |. Formally, we obtain this limit by taking 
Ad ^ 00 in eq. (10). Expanding 7± in powers of A±/Ad , 
7± 1 — A±/Ad, and eq. (10) yields: 



1 



< 



(11) 



The crossovers from attraction to repulsion are plotted in 
fig. 1 for several salt concentrations using eq. (10) and 
the no-salt limit, eq. (11). The second, well-known limit 
of high salt concentrations. Ad A±, was studied by 
Parsegian and Gingell [14] in the 1970s by linearizing the 
PB equation. They derived the pressure in the linear DH 
regime. In this case, 7± can be approximated by 7± ~ 
Ad/2A± in eq. (10), and the attraction condition is 



-d/Ar 



< 



< e' 



d/Xo 



(12) 



which reproduces the result of Parsegian and Gingell. In 
fig. 2 the line of zero pressure that separates attractive 
from repulsive interactions is plotted in the (a_/(7_|_, d/ Ad) 
plane for three salt concentrations. The plot shows the 
high salt limit of Parsegian and Gingell [from eq. (12)], 
as well as an intermediate amount of salt (A+/Ad = 0.3), 
and a relatively low amount of salt (A+/Ad = 0.05) from 
eq. (10). 
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Fig. 1: Regions of attraction and repulsion. cr+ and cr_ are 
the surface charge densities, and ad = l/2Trl bcI. The dashed 
and dash-dot lines represent the exact expression [eq. (10)] for 
crossover from attractive to repulsive interactions with d/ Ad = 
1 and 2, respectively. The solid line shows the no-salt limit 
[eq. (11)]. 




Fig. 2: The zero-pressure line 'n{a-/a+,d) = in the 
{a-/a+, d/Xo) plane. The dashed line represents the exact ex- 
pression [eq. (10)] for crossover with A+/Ad = 0.3, correspond- 
ing to Cb — 10^'' M for typical mica effective net surface charge, 
a+ ~ e/4000A^. On the dash-dot hue A+/Ad = 0.05, where 
for the same ~ 10~* M the surface charge is a+ ~ e/670 A^. 
The crossover in the DH limit [eq. (12)] is shown in solid line. 

We now compare the two figures and comment on the 
role of salt. In fig. 1 it is apparent that increasing the salt 
concentration enlarges the attractive region at the expense 
of the repulsive one. Another remark is that for |cr± | > ad 
(large d and/or strongly charged plates), the pressure is 
always negative (attraction) with no dependence on other 
system parameters like a-/a+ and Ad- This is related to 
the asymptotic behavior of the no-salt crossover curves. ^ 

^Another observation is that in fig. 2 it looks like the predictive 
behavior is opposite as compared to fig. 1 because for larger amounts 
of salt the repulsive region is enlarged. However, tfie two figures are 
in total accord because in fig. 2 the x-axis is also scaled by Ad . 



Gouy-Chapman pressure. — For the exact anti- 
symmetric case (cr+ =— (j_) and in the regime of low salt 
and strongly charged plates it was shown [15] that an ap- 
proximate scaling relation between the pressure 11 and 
d is: A/cf,/|n| ln(|n|/cf,) ~ d/Xn- The range of valid- 
ity of this scaling relation is the diagonal ray in fig. 3, 
(j+/(Jd = Ic-j/cTd ^ 1, and in addition d <C Ad- This 
result can be extended to the more general asymmetric 
region bounded by the hashed box in fig. 3: a+/ad 3> 1 
and \(j-\/ad ^ 1- From the integral in eq. (5) the relation 
between the pressure and d is deduced to be 




Fig. 3: The region of validity of the scaling relation is denoted 
by the hashed box in the {a+/ad, o-loa) plane. The dashed 
line denotes the anti-symmetric case. 

We note that in the entire range of validity of this scaling 
expression, the right hand side of eq. (13) varies between 
djXn and 3d/ Ad, meaning that it is roughly described by 
the anti-symmetric result mentioned above [15]. Thus, as 
long as the plates are strongly charged, n does not depend 
on the surface charge densities. This relation can be solved 
iteratively for 11 (rf) and the first iteration yields: 

n--^ln'W8AD) - (14) 

Counter-ion release. — In the limit of infinite sep- 
aration between the plates, an appropriate concentration 

of mobile cations and anions accumulate in the vicinity of 
each plate in order to neutralize the surface charge. When 
the plates are brought closer these two, oppositely charged 
clouds of mobile ions begin to overlap. Pairs of negative 
and positive counterions can thus neutralize each other 
and escape to the reservoir, where they gain entropy with 
no cost of electrostatic energy. This phenomenon of coun- 
terion release is the physical origin of the attractive forces 
between oppositely charged plates [15,16]. The parameter 
that characterizes the release of the counterions is defined 
as the excess charge per unit volume integrated over the 
entire separation between the surfaces:in eq. (10), 

r]= \ d0[c+(0)+c_(z)-2cfc]. (15) 

J-djl 
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where c±{z) are the number densities of the cations and 
anions. Safran [15] considered the relation between the 
pressure and the fraction of ions released into the reser- 
voir for the exact antisymmetric case, (cr+=— o-_), in the 
regime of low salt. Here we present a generalization of 
this result to the asymmetric case, where we focus on 
the same low salt, high surface charge regime but with 
Aa = a+ — \a-\ which is small compared to a± . Under 
these assumptions the integral in eq. (15) can be evaluated 
and gives: 



sA^ini^ - Act, 



(16) 



In the limit of 11 —> we find 77 ~ ct+ + |cr_ | ; this represents 
the largest possible value of excess counterions (beyond 
the bulk value ab = 2dcb). In this limit, corresponding to 
large separations [eq. (14)] each of the surfaces is neutral- 
ized by its own cloud of counterions. As the separation 
between the two surfaces decreases, the pressure reaches 
its maximal attractive value [II] 1/2ttIb^-, and the 
counterion excess approaches its minimal value: r] ~ Act. 
In this situation, not all the counterions are forced to re- 
main in the gap between the plates. The two oppositely 
charged plates screen each other except for an excess of 
surface charge Act that is compensated by the remaining 
counterions. 

The largest attractive pressure occurs for a distance d* 
that can be estimated as d* — 2A_ In 2Ad /A_ . Note the 
difference in this situation between the asymmetric and 
the exact anti-symmetric case for which the pressure sat- 
urates a,t d < d*. 

Scaling regimes. — In the general asymmetric case 
we must consider the charge asymmetry ratio, A-|-/A_ (or 
]cT_|/cT+) as an additional parameter. We investigate sev- 
eral scaling regimes in the three-dimensional parameter 
space: (A_|_/Ad, A+/A_ , d/Xu)- The pressure, as a func- 
tion of the asymmetry ratio, can be treated by consid- 
ering one of the following two limits: i) in the limit of 
\a-\/a+ <C 1, the negative plate can be taken as a neutral 
one, implying a repulsion between the plates, as demon- 
strated by eq. (10). In this limit, there is a mathematical 
correspondence between the asymmetric problem and the 
symmetric one (ct_|_ = ct_) with about twice the surface 
separation, d ^ 2d . The symmetric configuration satis- 
fies the boundary conditions corresponding to the case of 
one neutral and one charged surface by the vanishing of 
the electric field at the mid-plane. In this symmetric-like 
limit the pressure scales like the pressure in the symmetric 
CdbSQ-j clS discussed in detail at ref. [8]. 

ii) On the other hand, when the surface charge densities 
are nearly equal (and opposite), i.e. Act <C [cil , the pres- 
sure is attractive in a wide range of separations, and the 
formulae for the pressure are similar to those for the exact, 
antisymmetric case. Here, for simplicity, we present the 
regimes of the exact antisymmetric case with two oppo- 
sitely charged plates, ±ct. The different regimes are shown 
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Fig. 4: A schematic view of the various limits of the PB 
equation for two oppositely charged plates (the anti-symmetric 
case). The four regimes discussed in the text are: Debye-Hiickel 
(DH), Intermediate (Int.), Partial and Complete Counter- 
Ion Release (P-CIR and C-CIR). They are separated by four 
crossover lines: Int.^DH at Ago ~ Ad and Ad < d, DH<->C- 
CIR at \u = d and Ad < Ago, C-CIR^P-CIR at d = Ago 
and Ad > d, P-CIR<-^Int. at Ad = d and Ago < Ad. 



in the (Aqc/Ad, d/Xu) plane [fig. 4], and discussed below: 
(a) Debye-Hiickel. The limit of Ad «C Agc corresponds to 
low potentials (</> <C 1) for which the PB equation can be 
linearized. We obtain an attractive pressure (H < 0) as 
expected for two oppositely charged plates: 



H ~ 



(17) 



This pressure expression decays exponentially with dis- 
tance, (b) Complete Counter-ion Release. For small sepa- 
rations, <C Ad and d ^ Aqc , the charge neutrality of the 
system is maintained by the surfaces, and all the cations 
and anions are released to the reservoir. This yields a 
direct electrostatic interaction of two capacitor plates of 
charge ±cr in a dielectric medium: 



H; 



1 



2ttIbXq(^ 



(18) 



(c) Partial Counterion Release. This regime is defined by 
Agc ^ d <^ Ad, where the plates are strongly charged 
and the salt concentration is low. Equation. (14) can be 
used to predict the pressure as a function of the separa- 
tion, d. It is interesting to note that in this regime, the 
Gouy-Chapman pressure shows a long-range, power law 
dependence on d. It is also independent of the value of the 
surface charge, similarly to the result of Gouy-Chapman 
regime in the symmetric case, (d) Intermediate. When 
the plates are strongly charged, Xqc ^ Ad, the PB equa- 
tion cannot be linearized. However, if the separation is 
large, d ^ Xu, the surfaces are weakly interacting and 
can be treated as two separated plates. As a result, the 
electric field at the midplane is given by summing these 
two contributions. Under these assumptions, the electric 
field can be approximated by that of a single plate, and 
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the pressure is 

n^--^e-'^A-. (19) 

In the DH and intermediate regimes the pressure de- 
cays exponentially with the scaled separation, d/Ao due 
to screening effect of the salt ions. In the symmetric-like 
limit, the characteristic decay length is half as small. In 
the Gouy-Chapman regime the pressure is independent of 
the surface charge density, both in the antisymmetric and 
symmetric-like limits. This reveals a special property of 
the electrostatic interaction at medium-range separations: 
the length scales related to the surface charges, X± and 
Ad, have no effect on the pressure (beside a small logarith- 
mic correction as in eq. (14)). At small separations, the 
boundary conditions have the largest effect on the pres- 
sure between the plates. In the symmetric- like limit, the 
confinement of ions in between the plates results in a diver- 
gence of the; pressure, whereas in the anti-symmetric like 
limit, the pressure saturates due to the complete release 
of all the counterions at small separations. 

The similarity between the symmetric and anti- 
symmetric cases can be seen at large and intermediate 
separations, while the distinction is evident at small sep- 
arations. It is important to note that there is one major 
difference between these two limits: the symmetric- limit 
is purely repulsive, while the anti-symmetric one is always 
attractive. 

Discussion. — The theoretical investigation in this 
work suggests that the asymmetric Poisson-Boltzmann 
model predicts several interesting physical regimes in the 
interaction of dissimilar charged bodies immersed in elec- 
trolyte solutions. In particular, the pressure dependence 
on the separation Tl{d) can be evaluated from both the 
scaling relations and the numerical solutions. These re- 
sults can be tested in various experiments measuring forces 
between charged objects. One of the interesting results to 
be tested experimentally is the prediction for the crossover 
between attractive and repulsive interactions, 11 = 0, at 
high and low salt conditions. Others may include the 
crossover of the interactions from exponential decay to 
power law due to the onset of counterion release. 

We wish to point out two features of the behavior at 
small d. The first is that in our model a diverging repulsive 
pressure is obtained in the limit of vanishing d because of 
the assumption of fixed surface charge. However, effects 
such as charge regulation [17] and lipid demixing [18-20] 
can modify this assumption and may lead to to an overall 
attractive pressure for any d. 

The second is that attractive van der Waals (vdW) in- 
teractions always prevail at small separations, d < 2nm 
(in addition to the electrostatic interactions considered 
here). In experiments for small enough d, this vdW in- 
teraction is stronger than the electrostatic one and will 
overcome its repulsion. Thus, the only way to observe the 
pure electrostatic crossover between attractive and repul- 



sive interactions is to work in a set-up where this crossover 
occurs for a d range beyond the influence of vdW interac- 
tions. 
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